Weak antiferromagnetism and dimer order in quantum systems of coupled tetrahedra 
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We analyze the phases of an S=l/2 spin model on a lattice of coupled tetrahedra. The presence 
of both Heisenberg and antisymmetric, Dzyaloshinsky-Moriya interactions can lead to two types of 
symmetry-broken states: non-magnetic dimer order and, unexpectedly, exotic 4 sub-lattice weak 
antiferromagnetic order - a state with a generically small ordered moment and non-zero chirality. 
External magnetic field also induces weak antiferromagnetism co-existing with strong dimer corre- 
lations in the ground state. These states are formed as a result of broken Ising symmetries and 
exhibit a number of unusual properties. 
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I. INTRODUCTION 

Some of the most challenging and exciting problems 
in modern solid state physics are related to the na- 
ture of symmetry-broken states and the competition be- 
tween different types of order in insulating and doped 
antiferromagnetsi^ There exist numerous materials, 
ranging from the high-temperature superconductors to 
molecular-based magnets, that provide continuous source 
of inspiration for this research. 

In this paper we study the types of order that can oc- 
cur in quantum spin systems on lattices formed by cou- 
pled tetrahedra. Perhaps the most well-known (and still 
not fully understood) model of this type is the 3D py- 
rochlore lattice, composed of corner-sharing tetrahedra^ 
The purpose of this work however is to look at a class 
of models where the tetrahedra are connected weakly 
and in a more regular fashion, almost in a 2D square 
lattice- like arrangement. There are two main reasons 
for this. First, such models allow for more reliable and 
complete theoretical treatment. Second, we have had 
in mind potential applications to the 5=1/2 material 
Cu2Te2 05Br2* which is representative of such a geome- 
try. Rather unusually, in this material it has been ob- 
served that low-energy singlet excitations (measured in 
Raman spectroscopy) coexist with some kind of (pos- 
sibly weak) magnetic order, the origin of which is still 
controversial^ 

We suggest that to understand the properties of the 
above and similar systems, it is important to take into ac- 
count antisymmetric, Dzyaloshinsky-Moriya (DM) spin- 
spin interactions, which are expected to be present in 
a tetrahedron on symmetry grounds. Our main re- 
sult, which is quite general and model dependent only 
in the details, is that under certain conditions the low- 
energy singlet dynamics can coexist with weak antiferro- 
magnetism, induced by the DM interactions. It is well 
known that weak magnetic moments can appear near a 
magnetic-paramagnetic transition boundary, or in dimer 
(gapped) systems in the presence of both DM interac- 



tions and external magnetic field £ In the present work 
we describe a novel mechanism for weak antiferromag- 
netic order, which is induced by the DM interactions even 
without external fields. We show below that such an ex- 
otic possibility exists in tetrahedral systems due to the 
degeneracy of the ground state on a single tetrahedron. 
The typical excitation signatures are different from those 
of conventional spin waves and various nontrivial effects 
in an external magnetic field, such as magnetic field in- 
duced order, are also present. 

We start with the following spin-1/2 Hamiltonian, 
which involves both Heisenberg and antisymmetric ex- 
changes: 



W = ^JyS i .S j + ^D iJ .(S i xS j ), 



(1) 
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where the couplings Jij are distributed as shown in 
Fig. 1(b). The tetrahedra are represented as plaquettes 
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FIG. 1: (a) A tetrahedron with DM vectors represented by 
arrows (three of the six shown), (b) Two-dimensional lattice 
of coupled tetrahedra. 

in Fig. 1(a), and are assumed to be coupled weakly, i.e. 
_ <C J, where J is the exchange inside one tetra- 
hedron. All couplings are antiferromagnetic. We ana- 
lyze the model in two dimensions but the inclusion of 
three-dimensional couplings is straightforward and does 
not affect our main results. For simplicity of analysis 



2 



and notation we set for now j 2 = 0, and the effect of 
j 2 on our results will be discussed later. The DM cou- 
plings have relativistic origin^ and are expected to be 
weak, usually at most several percent of the Heisenberg 
exchange, |Dy| <C J. The distribution of the couplings 
Di j is very lattice-specific, as the presence of the DM in- 
teraction on a particular bond is determined by the sym- 
metry of the environment. We will make the assumption 
that the Di.j's are present only on the tetrahedra (where 
the Heisenberg exchange is also dominant), and are non- 
zero on every bond within one tetrahedron. The distri- 
bution of the six vectors Dij, written with respect to 
the principal axes of the cube in Fig. 1(a), and consistent 
with the tetrahedral group is: Di3 = ^=(— 1, 1,0), D 2 4 = 

(-1,-1,0), D 43 = ^(0,-1,1), D 12 = J|( ,-l,-l), 



D r(i>o, 



-1). Here D is the 



014 = ^(1,0,1), D 23 - 2 , 
magnitude of the DM vectors. 

The rest of the paper is organized as follows. In Sec- 
tion II we derive the effective Hamiltonian in the weak- 
coupling limit. The quantum phases of the model are 
then analyzed in Sections III and IV. Section V contains 
our conclusions. 



II. EFFECTIVE HAMILTONIAN IN THE 
PRESENCE OF DM INTERACTIONS 

First, we describe the states of one tetrahedron, 
Fig. 1(a). For D = the ground state is twofold de- 
generate (with energy Eq = — 3J/2) and is in the sin- 
glet sector S tot = 0. The two ground states are: |si) = 
^{[l,2][3,4] + [2,3][4,l]},|,s 2 ) ={[1,2][3,4]-[2,3][4,1]}, 

where [k, I] denotes a singlet formed by the nearest- 
neighbor spins k and I. The DM interactions break the 
continuous spin rotational invariance, leading to an ad- 
mixture of triplets to the ground stated Introducing the 
three excited (energy E% = — J/2) S = 1 triplet states 
Pn, In, tp, fx = x,y, z, in the notation of Ref. llOL we obtain 
the two new ground states, in the limit D/J <C 1: 

m = \si) + l\ Px ) - \ P y) + \q x ) + \qy)\ (2) 

I*) = 1*2) + 7^=J [\Px) + \Pv) + !<?*> - 1%)] + *-jl**>- 

These states remain degenerate with E$ M = — 3J/2 — 
3D 2 /(2J), and their wave-functions transform according 
to the e g irreducible representation of the tetrahedral 
point group. The diagonal matrix elements of on-site 
spins between the modified ground states @ are zero, 
whereas the off-diagonal matrix elements acquire finite 
imaginary values ( 1 3/|»S , Q i n 

The low-energy sector of the Hilbert space has dimen- 
sionality 2 tet , where N tet is the number of spin tetra- 
hedra. The remaining degeneracy is lifted by the inter- 
tetrahedral couplings, and we proceed to analyze the na- 
ture of the new ground state. For this purpose we in- 
troduce a pseudo-spin T = 1/2 representation, so that 



T z = 1/2 corresponds to |$) and T z — —1/2 corresponds 
to 1^). We then obtain the effective Hamiltonian in the 
ground state subspace: 



W eff = - 



<iJ> 



(3) 



Now the site indexes i, j refer to the positions of the tetra- 
hedra, and the summation is over nearest neighbors on a 
square lattice. To lowest order in ji,D the couplings in 
the different pseudospin directions are explicitly: 

7? AD 2 n 2 

n x = n y = j x — , n z = -f, h = -ci x . (4) 

The mixed coupling depends on the bond direction 
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FIG. 2: Phase diagram of the pseudospin Hamiltonian at zero 
temperature. 

The effective Hamiltonian describes an anisotropic fer- 
romagnet in an effective magnetic field h applied along 
the "hard" z-axis. If the field h does not exceed a critical 
value h c (which we find to be the physical case, corre- 
sponding to (@J), the system breaks the Ising symmetry 
along the soft x or y axes. The selected direction in 
the xy-plane depends on Q x /Q y . The first-order tran- 
sition line Q x — Qy, on mean-field level, separates re- 
gions with (T Xy i) ^ 0, (Ty.i) = (for Vl x > Q y ), and 
(T yA ) ^ 0,(T x 'i) = (for'f^ > Q x ). The two regions 
are shown in Fig. 2 as a function of the microscopic pa- 
rameters. Since h ^ 0, throughout the phase diagram 
(T 2j i) ^ 0. We now discuss in more detail the nature of 
symmetry breaking in the two phases. 



III. DIMER PHASE AND MAGNETIC FIELD 
INDUCED ORDER 

The phase with (T x \) ^ corresponds to dimer or- 
der. The existence of such a phase for the case D = 
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FIG. 3: (a) Dimer state, {T Xi i) 7^ 0. (b) Magnetic moments 
in an external magnetic field along the 1-3 bond in the dimer 
phase from (a) (tilt of spins in the field direction not shown), 
(c) Spin arrangement in the phase {Ty,i) 7^ 0, at zero field. 
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was pointed out in Ref. [llj. Physically, this means that 
the ground state of the system is the linear combina- 
tion ^|$) ± k\ty, on every site, which is equivalent to 
(T xA ) = ±V3/4, (T z ,i) = 1/4 (th ese values correspond to 
Q). The ground state is a real combination of |<fr) and 
|\&), and therefore the expectation values of local spins 
vanish, and only dimer order is present in this case, as 
shown in Fig. 3(a). Since the £l x coupling in (J3J) has 
ferromagnetic sign, the dimer pattern repeats itself on 
all tetrahedra. Magneto-elastic couplings could also con- 
tribute to the dimerization tendency as discussed for the 
pyrochlore case^ but in our scenario dimerization occurs 
spontaneously and is due to purely Heisenberg exchanges. 

The dimer ordered state breaks the discrete, Ising sym- 
metry T Xi i — > — T Xi i of ©, and thus the excitation spec- 
trum w(k) has a gap A = u>(k = 0) = 1.430a; (for 
D = 0). At finite temperature T, the ordered state ex- 
ists below the critical temperature T c — Q.92H, X . We have 
used the linear spin wave expansion and the mean-field 
equations to calculate these quantities. 

Due to the presence of triplets in the ground state , 
a (real) external uniform magnetic field — <?/ibH-S* ' also 
leads to non-trivial effects. On a single tetrahedron (ji — 
0) the field produces a finite magnetization in the field 
direction, shown in the inset of Fig. 4. For small fields 
H/J <C 1, the contribution is linear and proportional to: 
(D/J) 2 (H/J). We find that the magnetization is only 
weakly dependent on the field direction. For j\ 7^ 0, the 
corrections to the effective Hamiltonian Q have to be 
derived, and we find them to be strongly dependent on 
the field direction. For example for a weak field in the 
x-y direction H = ^=(1, 1,0), Eq. |3J) changes to: 



\*)h = \*)-^l\ qz )+V2(\t y )-\t x ))}. 



(5) 



The energy difference between these states leads to a 
modification of the "effective field" h -► h- D 2 H 2 / (2J 3 ) 
in Eq.(|3J). Even though several additional terms are gen- 
erated in H c ff, for j%, D, H -C J the above effect is domi- 



FIG. 4: Dependence of the gap (dashed lines) and T c (solid 
lines) on magnetic field (applied as in Fig. 3(b)) at the points 
A and B of Fig.2. Two upper curves shifted by 0.05 for clar- 
ity. Inset: Magnetic moment {Sl ot } in magnetic field in the 
z direction for a single tetrahedron (results identical in x,y 
directions). For D = 0, (Sf 4 ) = 0,H z /J < 1. 



nant. We have recalculated the gap A(H) and the critical 
temperature T C (H) in the presence of the field and the 
results are summarized in Fig. 4, where we have plotted 
these quantities relative to their zero field values. The 
gap generally shows a stronger field dependence. We 
note that the increase as a function of field is not uni- 
versal, and in fact if a field is applied in the z direction 
H = (0,0, if), then A(H) and T C (H) would decrease at 
almost the same rate. 

An applied magnetic field also generates staggered 
antiferromagnetic moments in the plane perpendicular 
to the field. The corresponding pattern is shown in 
Fig. 3(b), where (S 3 ) = -<S 4 ), (Si) = -<S 2 ), and the 
spins on sites 2 and 4 point along the diagonals of the 
cube (towards and out of the cube's center, respectively). 
The value of the magnetic moments in this pattern, re- 
peated on all tetrahedra, is 

\{S n )\ = ^-(T x , i ) = ^^, n= 1,2,3,4. (6) 

Field induced antiferromagnetic ordering is a general fea- 
ture of singlet systems with DM interactions. 



IV. BROKEN TIME-REVERSAL SYMMETRY 
STATE WITH WEAK ANTIFERROMAGNETIC 
ORDER 

The phase with (T y j) 7^ in Fig. 2 corresponds to a 
ground state which is a complex linear combination on 
a single tetrahedron a|$) ± where a, j3 are real 

coefficients. This combination is ferromagnetically re- 
peated on every tetrahedron. Depending on the values 
of the microscopic parameters, V3/4 < 1(1^)1 < 1/2. 
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The resulting state resembles to a large extent a spin- 
liquid state with a broken time-reversal symmetry pro- 
posed some time agc«i£ and characterized by finite scalar 
chirality \ = (Si • (S2 x S3)) 7^ 0. In our case such a state 
appears, however, due to the anisotropic DM interaction 
in a system without spin-rotational invariance, and we 
have x ~ (Ty,i)- Therefore, a broken time-reversal sym- 
metry immediately induces finite antiferromagnetic mo- 
ments with magnitude: 



|(S„)| = -V2(T y ,). (7) 



The spins are at an angle <j> = 109.47° with respects to 
each other, and form a four-sublattice antiferromagnetic 
structure, as shown in Fig. 3(c). On each tetrahedron 
the total magnetic moment ^2 n (S n ) = 0. 

Similarly to the dimer phase, the broken Ising symme- 
try leads to an excitation gap and a finite Ising transi- 
tion temperature, whose values, for dominant DM inter- 
actions are respectively A ~ Q y and T c ~ Q y . Thus, 
we find a quite unusual situation: the specific heat ex- 
hibits a large anomaly at the transition temperature, 
where a macroscopic part of the low-temperature entropy 
JVtet In 2 freezes out. However, if looked at in neutron ex- 
periments such a transition is characterized by develop- 
ment of small antiferromagnetic moments below T c . This 
exotic behavior is due to the fact that the transition it- 
self is driven by low-energy singlet degrees of freedom. In 
applied magnetic field H = 1) 0), the gap and T c in 

the broken time-reversal symmetry state scale as shown 
in Fig. 4, upper set of curves. We find that generally 
the field dependence is stronger than in the dimer phase, 
which is mainly due to the presence of the DM scale D 
in the two observables without a field. 

We note that Q y depends linearly on ji while J7 X is 
quadratic (0J. Thus, even though the T y (magnetic) and 
the T x (dimer) orders compete with each other, at suffi- 
ciently small ji the T y order is expected to dominate (al- 
though only in asymptotic sense for the present model). 
The tendency towards T y order is further enhanced by 
the presence of 7^ ( see Fig. 1(b)). This coupling frus- 
trates the T x component while strengthening T y , which 
results in a shift of the phase boundary in Fig. 2 down- 
wards. If we also consider a situation when ji is the 
dominant inter-tetrahedral exchange, J > J2 > ji, and 
set ji = 0, then we find that T x couplings are gener- 
ated only in fourth order, Q x ~ j%/ J 3 - This again means 
that the T y order can occur under much easier conditions. 
Nonetheless, only on lattices that have the tendency to 
produce true spin liquid ground states do we expect the 
T y order to occur spontaneously Good candidates ap- 
pear to be the pyrochlore lattice^ as well as the kagome 
latticed DM induced ordering has been discussed in the 
latter case from quasiclassical perspective^ 



V. DISCUSSION AND CONCLUSIONS 

Our calculations have been performed for ideal tetra- 
hedra while weak distortions (such as asymmetry in the 
Heisenberg exchanges, or deviations of the DM vectors 
from the tetrahedral symmetry) would introduce a small 
gap on a single tetrahedron. However, as long as j\,j2 
are sufficiently large, relative to this gap, the basic char- 
acteristics of the phases described above will remain un- 
changed. The DM order may deviate from Fig. 3(c), 
while preserving the property XL(Sn) = 0. 

The physics described in the present work occurs on 
the small energy scales f2 x , Cl v <C J, while upon increas- 
ing the inter-tetrahedral couplings to ~ J, a tran- 
sition to a more conventional Neel phase is expected to 
take place^i^ The exotic phases we have found are cer- 
tainly stable as long as the system is sufficiently far from 
the Neel order, i.e. the triplet gap is non-zero and all 
the relevant dynamics is governed by the singlet sector 
(or the mixed states Eq. (2) in the presence of DM in- 
teractions), and therefore can be described by the effec- 
tive Hamiltonian Eq. (3). We have made sure that this 
condition is fulfilled by using the weak-coupling analy- 
sis, although our results are expected to be qualitatively 
valid also for moderately large ji, J2 subject, as already 
mentioned, to the condition that the triplet gap is sizable 
enough. Thus, on general grounds we expect a transition 
as a function of j\ or j% between phases with exotic (weak 
magnetic) order, discussed in this work, and the Neel or- 
dered state, although the exact determination of the tran- 
sition boundary is a difficult problem and has not been 
the subject of this work. Instead, we have concentrated 
on the "universal" properties of the DM induced exotic 
phases which are clearly present in the weak-coupling 
limit and accurately described by a Hamiltonian of the 
form Eq. (3). 

The states we have found at low energies have quite 
unusual characteristics, which we will now summarize, 
(i) Due to the presence of the DM interactions the spin- 
rotational invariance is broken and the ground state de- 
velops a small magnetic moment in a field, with the uni- 
form magnetization showing a characteristic linear be- 
havior in weak magnetic fields (Fig. 4(inset)). (ii) Ei- 
ther dimerized or 4-sublattice, weakly antiferromagnetic 
states emerge as the ground state of the system. Both 
exhibit gaps in the spectrum since only Ising symme- 
tries are broken, (iii) The dominant components in the 
wave functions Eq.© are singlets, thus a large spectral 
weight S=0 transitions at low energy (equal to the gap 
~ Slz, fl y ) are present in the Raman spectrum, (iv) Both 
the gap and T c depend on magnetic field with a strong di- 
rectional dependence, (v) The dominant S=l excitation 
has a large gap of order J ^> £l x , fi y . 

We now comment on possible applications of our the- 
ory to Cu2Te20sBr2. The properties of this material are 
indeed consistent with signatures (i-iv) above, in particu- 
lar the compound exhibitsi*^ (a.) a weak magnetization 
in a field, (b.) a sharp zero-field, low-energy Raman peak 
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("singlet gap" A w 2AK), disappearing at T c w UK, 
(c.) a magnetic field dependence (increase) of both A 
and T c . Experimentally the most important unanswered 
question remains whether the ground state below T c is 
magnetically ordered or not, with recent NMR experi- 
ments pointing to some kind of magnetic ordertf This 
would suggest, in our scenario, that the ground state is of 
the type shown in Fig. 3(c). The crucial test for our ideas, 
and thus our prediction, would be the observation of a 
large triplet gap (larger than the "singlet gap" of 24K) in 
the neutron spectrum (point (v) above). The microscopic 
values of the exchanges are still rather controversial^ 
and consequently we have not attempted detailed fits 
but have used, phenomenologically, the model of Fig. 1. 
We emphasize that in our model a weak magnetic mo- 
ment is a generic feature of the (T y> i) ^ phase (and 
the dimer phase in magnetic field), and does not re- 
quire a fine-tuning of parameters to achieve proximity 



to a magnetic-paramagnetic boundary, as in the purely 
Heisenberg case£ 

Finally, on purely fundamental level, we have found an 
unconventional mechanism for weak antiferromagnetism. 
Our results are to be contrasted with the usual effect 
the DM interaction produces namely weak ferromag- 
netism in an otherwise antiferromagnetic system. In the 
present work we have shown how weak antiferromagnetic 
order can emerge from a singlet background. 
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